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ABSTRACT 



The nature of the interaction of atmospheric flow with a mountain 
range of finite length is investigated. An inviscid, adiabatic primi- 
tive equation model is used to simulate steady^ vertically-unsheared 
flow past a mountain range; various cases of mountain height and width 
are analyzed. A similar model for two-dimensional Boussinesq flow 
past a mountain on an f-plane is solved analytically following Merkine 
(1975) for comparison to the numerical solution. Results indicate that 
the atmospheric response is quite similar to that observed in the 
Boussinesq model. Ridging over the mountains and responses which damp 
in the vertical are observed in both models. Differences are noted in 
the magnitude of the response in the vicinity of the mountains and in 
the position of the downstream trough. In addition, the effect of 
finite difference model grid resolution is investigated by comparing 
simulations based on differing horizontal and vertical resolution. Re- 
sults indicate that for larger scale mountains, the effect of poorer 
horizontal resolution is significant; vertical resolution, however, is 
relatively unimportant. As the horizontal scale of the mountains is 
reduced, vertical resolution becomes increasingly more important so 
that the effects of both horizontal and vertical resolution are signi- 
ficant. 
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I . INTRODUCTION 



For many years, it has been clear that mountain barriers have impor- 
tant effects on cyclones and long waves as well as local effects such 
as mountain waves. The proper incorporation of topographic features and 
their effects in numerical models is, therefore, a prerequisite to 
accurate forecasts. 

Most numerical weather prediction models use a form of topography 
which is linearly smoothed according to grid size. Linear smoothing is 
a very serious restriction as may be seen in the hypothetical case posed 
by Charney (1967) . An infinitesimally narrow mountain ridge, oriented 
along a meridian and extending to the top of the atmosphere, would act 
as an absolute barrier to the flow. Smoothing, however, reduces the 
mountain height, and therefore its effect, to zero. 

An equally significant problem is that the dynamical mechanisms by 
which mountains influence large scale flow are still not fully under- 
stood. Earlier investigations, reviewed by Queney (1973), dealt with 
linearized flow over topography and were primarily concerned with 
smaller scale features. Recently, Merkine (1975) considered the inter- 
action of vertically-sheared, inviscid, steady flow with an infinitely 
long ridge; the formulation was based on approximating the horizontal 
momentum by the geostrophic momentum. Merkine and Kalnay-Rivas (1976) 
applied the same model to the three-dimensional flow over an isolated 
topographic feature. 

The main purpose of this study was to investigate the nature of the 
atmospheric response to mountains of varying widths and heights. The 
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primitive equation model developed by Monaco and Williams (1975) was 
used with a cyclic lateral boundary condition to simulate steady 
westerly flow past a mountain ridge of finite lengthy oriented perpen- 
dicular to the mean flow in the middle latitudes of the northern 
hemisphere. In addition, the ability of the finite difference model 
to resolve various features of the flow is considered. A simplified 
case with features quite similar to those of the model atmosphere was 
solved analytically in a manner similar to Merkine (1975). This solu- 
tion was then compared with numerical solutions obtained using various 
horizontal and vertical resolutions in the model. 
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II . DESCRIPTION OF MODEL AND INITIALIZATION PROCEDURES 



The two-level general circulation model used in this research has 
been described by Monaco and Williams (1975) . A brief summary of the 
model together with changes made during this research is presented here. 



A. MODEL DESCRIPTION 

The Monaco-Williams model is a modified version of the UCLA general 
circulation model described in detail by Arakawa and Mintz (1974) . The 
atmosphere is treated as adiabatic and frictionless with no sources or 
sinks. The equations are written in the vertical a-coordinate system 
where: 



a = 




( 2 . 1 ) 



7T E 



- P 



T 



where p is the pressure / p^ - is the constant pressure at the top of 

the model atmosphere, and p is the surface pressure. The vertical 

s 

boundaries of the model are the earth's surface and the 200 mb constant 
pressure level; hence, from (2.1), it follows that: 

0=1 at P = (2.2) 

0=0 at P = = 200 mb 

The boundary condition, O = ^ = 0, is applied at the upper and lower 
boundaries. The atmosphere in between is divided into layers of equal 
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Ac. The structure and the vertical distribution of variables is shown 



in figure 1 for the two-layer version of the model. The prognostic 
variables, which are horizontal velocity, temperature, and surface pre- 
sure, are governed by the primitive equations: 



Trr( u) + TrpC u) + -r— ( u) + u) - [ + (v u ) ] TTv 

9t mn 35 n 9n ni 9o mn mn 95 n 9n m 



7T,9$ 37T, tt „ 

+ -[^ + aa = — F- , 

n 95 95 mn 5 



(2.3a) 



9 ,iT , ^ 9 ,TTu . 9 ,7Tv . 9 ,7ra , ,f , 9 1 9 1 , , 

— ( — v) + — ( — v) + 3— ( — v) + 3— ( — v) + [ — + (v 33^ ^ 37 : “) 

3t mn 95 n 9n m 9a mn mn 95 n 9n m 



7r,9$ 9 tt, it „ 

+ — [t;;— + oa 3 —] = — F , 

m 3n 9ri mn r| 



(2.3b) 



3 /'n' 9 /'n'u 3 .TTv ,p ,K 9 .ira q, 

3 — ( — c T)+ 33 -( — c T)+ 3 — ( — c T) + (^— ) 3 — ( — c 0) = 

9t mn p 95 n p 9fi m p p 9a mn p 



,9 .TT , 



U 97T V 9tT, 77 

— ] + 

n 95 m 9ri mn 



(2.4) 



— (— ) 
9t^mn^ 



3 u 
n> 



d , V 
dT] m 



1 ^) = 


0 , 


(2.5) 


mn 


pa = 


RT , 


( 2 . 6 ) 


9<l> _ 


9a 


- TTa . 


(2.7) 



Here^ in ” ^ ^ ^ ^ complete listing of symbols for the above 

equations may be found at the beginning of this report. The model’s 
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diagnostic variables include the vertical a-velocity, o, and geopoten- 
tial, (p . The vertical Q-velocity is given by; 



a 



TT 

mn 



a 




( 2 . 8 ) 



o 



Geopotentials are determined in a manner which maintains the integral 
property of the vertically-integrated pressure gradient and conserves 
total energy. This requires combining the vertical integral of the 
hydrostatic equation (2.9a) with an alternate form of the hydrostatic 
equation written in terms of potential temperature (2.9b): 



The horizontal distribution of variables (figure 2) follows 
Arakawa's "Scheme C" which is designed to simulate accurately the geo- 
s trophic adjustment process. A detailed discussion of this scheme and 
its properties is given by Arakawa and Mintz (1974) . In this research, 
the latitudinal spacing of the TT points was fixed at 4®. In the 
longitudinal direction, two different grid resolutions were considered: 
coarse (AX = 4.5°) and fine (AX = 1.5°) . Cyclic continuity with wave 
number 8 was assumed in the ^-direction. To avoid the use of an ex- 
tremely short time step due to the convergence of meridians at the 
poles, the ^-component of the pressure gradient force and the divergence 
in the ^-direction are zonally smoothed as noted by Monaco and Williams 
(1975). The time integration proceeds in sequence of five time steps: 
one Matsuno forward/backward step followed by four leapfrog steps. As 



1 



1 




(2.9a) 



(2.9b) 
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7T 7T • 

noted by Williams (1976)^ the definitions of ( — ) and ( — O) at u and v 

mn mn 

points in the finite difference forms of the zonal and meridional momen- 
tiim equations were changed to the more complex forms described by 
Arakawa and Mintz (1974) . 

B. ANALYTIC INITIALIZATION 

Initial conditions were similar to those described by Monaco and 
Williams (1975) . Haurwitz (1940) and Neamtan (1946) obtained the stream 
function solution to the linearized barotropic vorticity equation. 
Phillips (1959) used this solution in the forcing function to obtain 
the geopotential from the non-linear balance equation. 

Following these results, an atmosphere in solid rotation (no vertical 
shear) was simulated by initializing the model with the following 
equations : 

2 

= - B a sin (p 

$' = i B a^(2Q + B) (2.10) 

1 dllJ „ , 

u = - — = Ba cos (p 

a d(p 

1 d'P ^ 

r = 0 

a cos (p dA 

where is the geopotential perturbation and B is a constant; the 

other terms are defined at the beginning of this report. By setting 
20 -1 

B = — sec , a zonal wind distribution which varies as the cosine from 
a 

the equator (where u = 20 m-sec ^) to the poles (where u = 0 m-sec ^) 
was initially specified. The initial temperature distribution was accord 
ing to the NACA standard atmosphere as defined in Haltiner and Martin 
(1957). 
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C. INTRODUCTION OF TOPOGRAPHY INTO THE MODEL 



An initial problem to be overcome was the adjustment of the model to 
the topography. Adjustment (or, in this case, steady state) was assumed 
if all surface pressure tendencies were of the order of 0.25 mb/hour or 
less and fluctuations in certain arbitrarily selected wind components 
were of the order of 0.1 m-sec ^/hour or less. A special scheme was used 
to introduce the topography into the model in order to minimize the time 
required to reach steady-state. The height of the lower boundary was 
initially set to zero everywhere. Heights at mountain grid points were 
then incremented during the integration as a function of time until the 
desired heights were reached. The equation used at these grid points is: 



$ = 
s 



K. • 2 ,7Tt 

M sin 

2t 



t < t 
- P 



( 2 . 11 ) 



M 



, t > t 
- P 



where is the height of the mountain at time, t; M is the desired 

mountain height; and t^ is the period of time over which the mountains 
are "built". A typical time series of the surface pressure tendency at 
a grid point in the vicinity of the mountains is shown in figure 3 for 
an integration of the two- layer fine resolution model with t^ = 16 hours. 
Note that the largest tendency (~15 mb/hour) is observed at mid-period; 
this corresponds to the maximum rate of change of the mountain height de- 
fined by (2.11). With one exception, it appears that steady-state has 
been achieved at 16 hours. One attempt to reduce t^ to 12 hours resulted 
in lengthening the adjustment time. In addition, it was necessary to 

increase t to 36 hours for the coarse resolution model. 

P 



17 



► n 




III. ANALYTIC INVESTIGATION 



Merkine (1975) considered a steady, Boussinesq flow with no vertical 
shear past a large mountain ridge in a rotating, stratified atmosphere, 
bounded above by a horizontal rigid lido He was able to analytically 
solve the potential vorticity equation with appropriate boundary condi- 
tions; he was also able to derive and solve an analytic equation for the 
V component of the wind. He performed this analysis on an f-plane with 
infinite boundaries in the upstream and downstream directions; in addi- 
tion, he assumed the flow to be balanced in the vertical and cross- 
topography directions (similar to Hoskins and Bretherton, 1972; and 
Hoskins, 1974) . If these same assumptions are made together with a 
lateral boundary condition of cyclic continuity, an analytic solution 
can be derived for a simplified case which is quite similar to the 
numerical integrations which are presented in this report. The two- 
dimensional potential vorticity equation (independent of y) for inviscid, 
adiabatic, Boussinesq flow on an f-plane may be written: 




9v 

9x 



+ f) 



3z 



3v 36 
^ 3x^ 



0 



(3.1) 



Consider the solution: 



. 3v ^ ^ 

3x 3z 3z 3x 



= constant 



(3.2) 



In order to determine the constant, equation (3.2) is evaluated away from 
the mountains where there is no horizontal or vertical shear to give: 

30s 

f — = constant (3.3) 

oZ 
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where 9 (z) represents the mean state potential temperature at level z. 

1 8$ 

using the geostrophic relation, v = — r and the hydrostatic relation, 

0 

o <3 V 

0 = 9 (z) + — - 77 - , equation (3.2) may be written: 
s g dz 



1 3^ „2 , 1 1 ,3^t >2 _ „ 

fTT'' ^ -2 TT rr 7T - -2 < 3Si-> ■“ 

3x f dx dz dz f 



(3.4) 



o 99 

2 Q s 

where N = 7 ^ 



0 3z 
o 

In order to non-dimensionalize this equation, the following scale 
factors (also used by Merkine (1975) ) are introduced: 



$ = 



2 2 
N H 



z = 



X = 



Hz* 

m 



X* 



(3.5) 



The non-dimensional form of (3.4) may then be written: 



d^4>' 






3x 



2 2 2 

*^ 3x* 3z*^ 



2 2 
+ S _ ( 3 ^ 



3z' 



3x * 3z 



r) = 0 , 



(3.6) 



or, dropping the prime notation: 



$ (1 + 0 )+0 -0 
XX zz zz xz 



= 0 , 



(3.7) 



where the subscripts x, z indicate differentiation with respect to a 

30 

particular variable - e.g., 0^ = ^ . This is similar to the equation 
derived by Merkine (1975) ; he further showed that this equation can be 
transformed into Laplace's equation under a non-linear transformation. 
First, equation (3.7) is differentiated with respect to z , and (p is 
replaced by through the following equation: 



ip = 0 + z - 1 

z 



(3.8) 
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Here, ijj is the non-dimensional stream function for flow in the x-z 



plane. In a steady-state, the non-dimensional potential temperature, 

0 + z , is constant along a streamline. The new equation in terms of 

z 

\p can be written: 

\p lb ^ -I- (lb ^ + 1) ip - 2ib lb lb =0. (3.9) 

^XX Z X ^zz ^x^z^xz 



Now a transformation of variables is made so that \p is the new verti- 
cal coordinate. With this transformation equation (3.9) becomes: 



z 

XX 






0 , 



(3.10) 



where z is now a dependent variable which represents the height of 
the isentropic surfaces. It is convenient to define the new variable 



U = z - \p -h 1 . (3.11) 

In terms of ^ the basic equation and boundary conditions are: 



XX ipip 



= 0 



y = 0 at \p = 0 







(3.12) 



1/2 

where f(S x) is a function specifying the height of the lower 
boundary. Since the numerical model invokes cyclic continuity as the 
lateral boundary condition, let 



f(S^^^x) = cos (S^^^x) 



(3.13) 



27T T 

, „l/2 R 2Triffl . ^ 

where S = — = — — ; in dimensional form (3.13) becomes 

Li lLj 
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f 

I 






27TX* 

COS ( — £ — ) t ^here the (*) denotes dimensional form. Setting y = 

1/2 

cos (S x)G(^) f it may be shown that (3.12) reduces to: 

- SG(\p) = 0 

G(0) = 0 , (3.14) 

G(-l)= I 

which has as a solution: 

1/2 1/2 

G(ip) = A sinh(S ip) + B cosh(S ^ Jp) . (3.15a) 



Equation (3.15a) is evaluated at the boundaries to determine the constants 



A and B with the following result: A = 

the solution to (3.12) may be written: 



h 

1/2 

H sinh S ^ 



B = 0 ? hence 



y 



h 

H 



sinh 

sinh(S^^^) 



/oV2 , 

cos(S x) 



(3.15b) 



Making use of the transformation rules and the equations of motion^ 
it can be shown that: 



3v _ 

d\p 9x 

, (3.16) 

9y _ _ 

9x dip 



and that the v component is determined by the system: 



V 

XX 



+ V 






0 



V 

V 






= 0 at 



^ H ^ 






1/2 



= 0 



f (S^/^x) 



at 



4 ; = -1 



(3.17) 
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this 



1/2 1/2 1/2 
By letting f ' (S ' x) = - sin(S x) and v = F(\p) sin(S x) , 



system may be solved in the same manner as above; hence, 

,1/2 



v(x,\p) = - 



h cosh(S 



H sinh S 



1/2 



sin(S^'^^x) . 



(3.18) 



Note that 



V = v(x,0) 
top 



H sinh S 



1/2 



sin(S^'^^x) , 



(3.19a) 



h . _l/2 , 

V, . . = v(x,-l)= :t7T sin(S x) , 



bottom 



H tanh S 



1/2 



(3.19b) 



are the solutions at the upper and lower boundaries. The fully dimen- 
sional forms of V at the lower and upper boundaries may be written: 



"b* = 



Nh . ,27 Tx*, 

sin (— ) , 



, ,2TTNH, 
■tanh(— ) 



(3.20) 



= 



Nh . ,27rx*. 

sin (— ] 



. , ,27TNH, 
sinh (— r ; — ) 



(3.21) 
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The asterisks are used here as a reminder that x and v are in their 



dimensional forms; represent values at the lower and upper 

boundaries respectively. Considering the lower boundary, v^^* may be 
written as the sum of an infinite series: 

Nh 



V, 



= - z 



^ , . 27TNH , 

n=l 



n . ,27Tx* 

sin (• 



^/n • 



(3.22) 



In a finite difference model, this summation must be truncated after 
wave number = 2 ^ (the maximum wave number resolvable by the grid) 
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The height of the mountain ridge, h^ , is specified in the following 
manner. In the dimensional form, , the height of the lower 

boundary may be written: 



z, * = h cos (— ) 

D L 

= h f(x*) , 



(3.23) 



where f(x*) is a functional representation of the shape of the topo- 
graphy varying from 0 at the base of the mountain to 1 at its top. 

In this investigation, only mountains symmetric (in the east-west 
direction) about the ridge line were considered (see figure 4) . The 
exact functional form of the mountain may be written: 



f (• 



nx*^ 



00 



1 

0 



-E 



n I X 
<£ 7T£ 



- eiT _< X _< eiT 
ETT _< |x| _< IT 



(3.24) 



(3.25) 



nx* 

f (-^ — ) ^Iso be represented by the following Fourier-cosine series: 

00 

, ,27Tnx*, 

f (-^~) = + y An cos (- ^^ — ) , 

e 

where: A = — , 

2 (1-cos (neiT) ) 



A = 
n 



2 2 
n 7T e 



(3.25a) 



The height of the mountain ridge (at x* = 0) may then be written: 





00 






h = h + h , 

o jLj n ' 


(3.26) 




n=l 




where: 


h = h A 


(3.26a) 




o o 




h = h A 
n n 


(3.26b) 



h^ does not contribute to the v component; therefore, (3.26b) is used 

in equation (3.22) to derive v. * . 

b 
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In order to analytically approximate the conditions present in the 
numerical model, the series approximations for v^* and were 

truncated after wave number n = • For example, a 30-point grid 

can resolve only 15 wave numbers; therefore, the series approximations 
are truncated after 15 terms. The v^* analytic profile for a 
9® X 1.5 km (width x height) mountain computed on a 30-point grid span- 
ning one cycle of wave number 8 is shown in figure 5 where the origin 
of the graph is placed at the mountain top. Since a 30-point grid on 
wave number 8 results in 1.5° resolution, the mountain is carried (by 
non-zero at five grid points. Superimposed is the same analytic 

case computed with twice the resolution; this results in the mountain 
being represented by 11 grid points. Note that the two profiles at the 
lower boundary are nearly coincident. 

In figure 6 is shown the v^^* analytic profiles for a 3° x 1.5 km 
mountain computed using a 30-point grid versus the same case using a 
60-point grid. In the former case, only one grid point represents the 
mountain, while three points were used in the latter case. Qualitatively 
the two profiles exhibit the same general flow characteristics; however, 
the fluctuations observed in the profile of the 30 grid point case 
differ significantly from the smooth profile of the 60 grid p6int case. 
When the mountain width was reduced to 1.5° on the 60-point grid (there- 
by allowing only one grid point to represent the mountain) , similar 
fluctuations were observed. These results seem to indicate that one grid 
point is an inadequate representation of a mountain with respect to the 
quantitative aspects of the flow at the lower boundary. 

In figure 7, the v^* analytic profiles for a 3° x 1.5 km mountain 
computed on a 30-point grid versus a 60-point grid are shown. The 
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similarity between the two profiles indicates that the poor resolution 
has little effect at the upper boundary. As noted above, a grid with 
1.5® resolution appears to be sufficient to accurately specify the flow 
response to a 9® wide mountain at both the upper and lower boundaries. 
Therefore, in order to examine the characteristics of the analytic re- 
sponse, consider figure 8 in which the analytic profiles of and 
are shown for a 9® x 1.5 km mountain computed on a grid with 1.5® reso- 
lution. Both profiles are antisymmetric about the mountain indicating 
ridging centered at the mountain top. A downstream trough, commonly 
associated with flow over mountains, is located at the eastern boundary. 
However, in this case, the feature is purely the result of the assumption 
of cyclic continuityo Had we assumed the lateral boundaries to be in- 
finite, we would have observed downstream a permanent deflection to the 
south (as noted by Merkine (1975)) because of the assumption of a constant 
Coriolis parameter^ Merkine (1975) also noted that a variable Coriolis 
force would generate a trough in the downstream region (results from the 
numerical model support this observation). The lower boundary response, 
as expected, is greater in magnitude than the upper boundary response. 

Note that the lower boundary maxima, however, are located closer to the 
ridge. This would seem to indicate that the smaller scale response is 
damped more rapidly above the lower boundary. At the ridge itself, both 
v^* and v^* are zero. It may be seen from (3.20) and (3.21) that varying 
the mountain height affects the amplitude of the response only. This 
can be seen from figure 9 in which the v^* profile is shown for a 9® x 
1.5 km mountain vs. a 9® x 3.0 km mountain. Increasing the height of the 
mountain by a factor of two causes a uniform increase (by a factor of two) 
in the magnitude of the response at both the upper and lower boundaries. 
This is not true in the interior where non-linear distortions occur. 
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IV. NUMERICAL MODEL RESULTS 



A number of integrations were performed using various forms of the 
model. In each of the experiments, a mountain range of uniform height 
oriented north-south was centered in the grido The model was initialized 
as earlier described, and integrated forward in time until a steady-state 
was achieved. In the following discussion again, fine resolution refers 
to 1.5° longitudinal grid spacing and coarse resolution refers to 4.5° 
grid spacing. 

Ao BROAD MOUNTAIN RANGE CASE 

A mid-latitude mountain range, 9° x 1.5 km (width x height), extend- 
ing from 34°N to 66°N was the first case considered. Results from the 
various integrations are as follows: 

1. Six Layer Fine Resolution Model Results 

The wind field from level six (O = .91667) is shown in figure 10. 
Significant features include ridging over the mountains and weak trough- 
ing downstream. Upstream, the flow is deflected northward by the 
mountains; downstream, the deflection is southward, becoming weakly north 
ward east of the trough. At a distance from the mountains, both up- 
stream and downstream, little deflection is apparent. Deflection of the 
flow is also observed one grid point (-445 km) to the north and south of 
the mountains; this distance is of the order of the Rossby radius of 
deformation. In figure 11 is shown the wind field at level one (a = 
.08333) . In general, the effect of the mountains on flow at this level 
appears to be minimal. Consider the v profile for level six (0 = .91667) 
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at 44®N shown in figure 12. This latitude was chosen because it crosses 
the mountain near the center and, more importantly, because it was re- 
presentative of the response at the other latitudes which cross the 
mountains. Superimposed on the v-profile is the analytic solution, f 
at a = 1.0 for the similar analytic case. There is significant similar- 
ity between the two curves. Both profiles indicate ridging at the 
mountain top with the only apparent upwind difference being one of magni- 
tude (the model showing a greater northward deflection). Dov/nstream, 
as in the analytic case, the southward maximum is in the same position, 
relative to the mountain top, as the northward maximum; however, its 
magnitude is about twice that of the northward maximum. Also, downstream 
the model profile indicates the trough position to be about 75° (longi- 
tude) east of the mountains; in contrast, the analytic solution shows 
the trough at the channel *s eastern boundary. This difference is to be 
expected because the Boussinesq approximation and the f-plane and rigid 
lid assumptions (present in the analytic case) are not made in the 
numerical model. In figure 13, the v profiles for levels two, four and 
six at 44°N are shown. Noting the positions of the northward and south- 
ward maxima relative to height, a systematic westward tilt in the verti- 
cal is observed; this tilt is not observed in the analytic solution. 

It is felt that this tilt could be the result of the 3-effect and other 
effects not included in the analytic solutiono As noted in the analytic 
case, the response should damp with height — this is especially apparent 
in the magnitude of the southward maximum. In figure 14 is shown the 
analytic solution at a = 0 versus the model profiles for a = .08333 and 
O = .25. Note that in both model profiles the magnitude of the response 
is similar to that of the analytic case, but the horizontal scale is 
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smaller. In addition, the systematic westward tilt with height is ob- 
served in the two model profiles. 

2. Two Layer Fine Resolution Model Results 

The same case was computed using the two layer model to examine 
the effect of vertical resolution on the numerical simulation of the 
flow. The a = *75 velocity field was compared with its counterpart from 
the six layer model and only minor differences (in the vicinity of the 
mountains) were observed. In figure 15, the v-profiles for a = .75 at 
44°N are shown for these two models. There is little difference between 
the two curves, both showing a ridge centered one grid point west of the 
mountains and both showing weak troughing in the same general area 
downstream. Some differences exist in the vicinity of the mountains, 
but, in general, the two layer model's response appears to be quite 
similar to that of the six layer model. In figure 16, the v-profiles at 
a = .75 and a = .25 from the two layer model are superimposed. As in 
the case of the six layer model, a westward tilt is again observed. 

3. Coarse Resolution Model Results 

In order to investigate the effect of horizontal grid resolution, 
the above integrations were repeated using the six-layer and two-layer 
coarse resolution models. The velocity fields at a = .25 and a = .75 
produced by these two models were quite similar. In figure 17 is shown 
the six-layer model field at a = .75. In general, it shows a deflection 
of the flow similar to the six-layer fine resolution model, but differ- 
ences were noted particularly in the vicinity of the mountains. In 
figure 18 are shown the v-profiles for the six-layer fine resolution model 
versus the six-layer coarse resolution model. Note that the magnitudes 
are about the same but that the scale of the response is less in the 
fine resolution model. This would appear to emphasize the importance of 
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horizontal grid resolution in the simulation of finer scale features 
of the flow over topography. 

Consider figure 19 in which the lower level (a = .75) v-profiles 
for the two-layer coarse resolution model are shown for three different 
mountain height cases (1.5, 2.4, and 3.0 km). The shapes of these 
three profiles are quite similar with magnitudes increasing with in- 
creasing mountain height. Clearly, the response is directly proportional 
to mountain height as it is in the analytic case at the boundaries. In 
addition, although the response at the mountain top is non-zero, it is 
approximately independent of mountain height. Hence, at least some of 
the gross aspects of the analytic solution are apparent in the coarse 
resolution simulations. 

B. NARROW MOUNTAIN RANGE CASE 

Next, a narrower mountain range (3° wide) extending from 34®N to 
SS'^N was considered; mountain height was again set at 1.5 km. Since 
the distance between grid points in the ^-direction is 1.5® in the fine 
resolution model, this allows only one grid point at each latitude to 
represent the mountains. 

1. Six Layer Fine Resolution Model Results 

The wind field from level six (O = .91667) is shown in figure 
20. As in the wider mountain case, there is ridging over the mountains 
and troughing downstream. However, the magnitude of the deflection in 
the vicinity of the mountains appears to be much greater in this case. 
Note that one grid point (1.5® longitude) west of the mountains, the 
flow is almost directly from the south the whole length of the mountain 
range. Significant deflection of the flow is observed one grid point 
(^445 km) north and south of the mountains as was true in the wide 
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mountain case. The effect of the mountains on the upper level flow at 
a = 0.08333 appeared to be minimal (as in the wide mountain case) . In 
figure 21 is shown the v-profile for level six (a = 0.91667) at 44°N 
superimposed on the analytic solution at a = 1.0 for the same case. The 
strong northward deflection observed in figure 20 is seen here at grid 
point 14 as a strong northward maximum exceeding, in magnitude, the 
analytic solution by a factor of four. While the downstream southward 
maximum is not coincident with its analytic counterpart, it has the same 
magnitude as the northward maximum. This symmetric behavior is predicted 
by the analytic solution and it differs from the result in the wide 
mountain case where the magnitude of the southward maximum was about 
twice that of the northward maximum. 

The position of the downstream trough is closer (by 3*^ longitude) 
to the mountains in this case. Note the similarity between the fluctua- 
tions observed in the analytic solution and those appearing in the up- 
stream v-profile. It would appear that these oscillations are a result 
of allowing only one grid point to represent the mountains interacting 
with a mean flow from the west. These oscillations are not present in 
the model’s downstream profile. In figure 12 are shown the v-profiles 
for levels two, four and six. Noting the position of the southward maxi- 
mum with respect to height, there again appears to be a systematic west- 
ward tilt with height; however, it is not as well defined here as in 
the wider mountain case (figure 15) . In figure 23 is shown the analytic 
solution at a = .08333. As in the wide mountain case, the amplitudes of 
the two responses are similar, but the horizontal scale of the model re- 
sponse is much smaller. 
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2 . Two Layer Fine Resolution Model Results 



An integration for the same case was performed using this ver- 
sion of the model. Comparing the resulting wind field at a = .75 with 
the six layer model results, significant differences were observed in 
the vicinity of the mountains (see figure 24) . In the six-layer model 
field, a well-defined, narrow trough broadening toward the south is ob- 
served immediately east of the mountains; however, there is no trough 
in this area in the wind field of the two layer model. Differences 
were noted in the magnitudes of the wind velocity downstream. Fur- 
ther upstream and downstream, the two models give similar results at 
this level. A comparison of the fields at a = .25 showed little differ- 
ence between the two models. In figure 25 are superimposed the v 
profiles from both models for a = .75 at 44®N. As can be seen, little 
differences exist away from the mountains; however, immediately upstream, 
the oscillations in the six-layer model profile appear to be of greater 
amplitude, and, downstream, the six- layer model profile shows the 
presence of a strong trough which is absent in the two-layer model pro- 
file. In figure 26, the v-profiles from both models at a = .25 (also 

at 44°N) are shown; the two are similar to the extent of having a prin- 

% 

cipal trough and ridge, but sizable phase differences are seen. These 
results indicate that the vertical resolution is important in this case. 

During the latter stages of this study, a relatively minor coding 
error was discovered in the computer program; this error had been 
present throughout the period of this research. As described by Monaco 
and Williams (1975) , a pressure-area-weighted form of the horizontal 

TTU 7TV 

Wind components ( — , ) is forecast by the model; u and v are computed 

mn mn 

TT 

after each time step by dividing out — . The error resulted in a slight 



overestimation of this term in that section of the coding responsible 
for removing this weighting from the v component after each leapfrog 
time step. A test case was recomputed using the corrected version of 
the two-layer coarse resolution model. Almost no change was observed 
in the wind fields for the Q = .25 and a = .75 levels. Although slight 
changes in the v-profiles were observed, the overall pattern remained 
the same. It was, therefore, concluded that the error in coding did 
not significantly affect the nature of the response depicted in these 
results . 
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V. CONCLUSIONS 



From the results of the numerical simulations, it is evident that, 
in general, the mountains interact with the atmospheric flow in such a 
way as to induce ridging over the mountains and troughing on the down- 
stream side at lower levels. With the exception of the troughing the 
numerical solution appeared, in general, to agree with the analytic 
result from the Boussinesq model. The presence of the downstream trough 
was due to the fact that the numerical model has a variable Coriolis 
parameter and a non-rigid upper boundary; in contrast, the analytic 
model considered Boussinesq flow on an f-plane bounded above by a non- 
rigid lid. 

Several other similarities were observed between the analytic 
response and the numerical response. As in the analytic case, smaller 
scale features appeared to damp rapidly with height. Varying the 
mountain height was, in general, observed to affect the amplitude of the 
numerical response in the same manner as in the analytic case. Reduction 
of the width of the mountains reduced the horizontal scale of the 
numerical response, but in contrast to the analytic solution, it also 
increased the magnitudes of the northward and southward maxima. It was 
initially felt that the strong northward maximum along the western side 
of the 3® X 1.5 km mountains might be related to air flowing around 
rather than over the mountains. A detailed examination, however, 
seemed to indicate that this was not the case. 

The numerical solution, in every case, displayed a systematic west- 
ward tilt with height which was not observed in the analytic solution. 
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The precise reason for this tilt was not determined; however, it appears 
to be the result of the 3*"effect, compressibility and non-rigid upper 
boundary in the niimerical model. 

The effect of vertical resolution appeared to be minimal in the 
case of the broader mountains; but significant effects were observed 
when the mountain width was reduced. With respect to horizontal resolu- 
tion, the coarse resolution model appeared to be able to reproduce some 
of the larger scale features observed in fine resolution solutions; 
however, differences were noted between the two models with respect to 
finer scale features of the flow. 

Early in this research, an unsuccessful attempt was made to simulate 
conditions in which air flows around rather than over the mountains. 

From these results, and those of Merkine (1975), it appears necessary 
to further reduce the horizontal scale of the mountains. These results 
also indicate that finer horizontal and vertical resolution would be 
necessary to accurately simulate the atmospheric response to this 
narrower mountain range o With the high price paid in computing time 
for increased resolution, this appears to be a definite argument for 
the use of a nested grid in treating atmospheric response to mountains. 

No problems were encountered in treating unsmoothed topography in 
the numerical model. In addition, the scheme used to introduce the 
mountains into the model behaved quite well. Additional attention 
should be given to the possibility of reducing the period of time neces- 
sary to incorporate the mountains. If this can be done, a similar 
scheme might be used to introduce smaller scale topographical features 
into operational numerical models. 
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Finally, this investigation considered a mean flow with no shear in 
the vertical* Future research should be directed at investigating the 
response of a vertically-sheared flow to mountains. Another area to 
be examined is the effect of condensation and the release of latent heat 
into the air flow. 
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Figure 1. The sigma (a) coordinate system as used in the two-layer 
version of the model. Also shown is the vertical distri- 
bution of the variables. 
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Figure 3. Time variation of the surface pressure tendency from a 
grid point in the vicinity of the mountains. 
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Figure 5„ 9° x 1,5 km mountain case. Analytic v profiles for a = 1.0 

at 44°N computed on a 30 point grid (dashed curve) vs. a 60 

point grid (solid curve) . Vertical axis is centered at the 
mountain top. 
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Figure 6. 3® x 1.5 km mountain case. Analytic v profiles for O = 1.0 

at 44®N computed on a 30 point grid (dashed curve) vs. a 60 
point grid (solid curve) . 
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3° X lo5 km mountain case. Analytic v profiles for 0 = 0.0 
at 44®N computed on a 30 point grid (dashed curve) vs. a 60 
point grid (solid curve) . 
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Figure 8. 9® x 1.5 km mountain case® Analytic v profiles for o = 1.0 

(solid curve) and o = 0.0 (dashed curve) computed on a 30 
point grid® 
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Figure 9. Analytic v profiles for a = 1.0 at 44 °N computed on a 

30 point grid for a 9° x 1.5 km movintain (dashed curve) 
vs. a 9° X 3,0 km mountain (solid curve). 
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Figure 10. 9° x 1.5 km mountain range case. Six-layer fine resolution 

model streamlines for o = .91667. Dashed lines depict out- 
line of mountain range and mountain ridge line. 
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Figure llo 9® x I® 5 km mountain range case* Six— layer fine resolution 
model streamlines for a = « 08333. 
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Figure 12. 9° x lo5 km mountain range caseo Analytic v profile for 

CJ = loO (dashed curve) vs. six-layer fine resolution 
model V profile for o = .91667 (both at 44®N) . 
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Figure 13. 9° x 1.5 km mountain range case. Six-layer fine-resolution 

model V profiles for a = .25 (dotted curve), O = .58333 
(dashed curve), a = .91667 (solid curve), all at 44°N. 
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Figure 14. 9° x 1.5 km mountain range case. Analytic v profile for 

a = 0.0 (dashed curve) vso six-layer fine-resolution 
model V profiles for a = o083333 (solid curve) and 0 = .25 
(dotted curve) , all at 44®No 
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Figure 15. 9° x 1.5 km mountain range case, v profiles for o = ,75 at 

44®N from six-layer fine resolution model (solid curve) and 
two-layer fine resolution model (dashed curve) . 
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Figure 16o 9° x 1.5 km mountain range case, v profiles for a = .25 
(dashed curve) and a = .75 (solid curve) at 44®N from 
two-layer fine resolution model. 
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9® X 1.5 km mountain range caseo Six-layer coarse 
resolution model streamlines for a = .75. 
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Figure 17. 




Figure 18. 9® x 1.5 km mountain range case^ v profiles for O = .91667 

at 44®N from six-layer fine resolution model (solid curve) 
and six-layer coarse resolution model (dashed curve) o 
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Figure 19. v profiles from two-layer coarse resolution model for 
a = o75 at 44®N for the following cases; 9° x 1.5 km 
mountain range (dotted curve ) f 9° x 2.4 km mountain range 
(solid curve), 9® x 3.0 km mountained range (dashed 
curve) o 
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Figure 20» 3° x 1»5 km mountain range case.. Six-layer fine resolution 

model streamlines for a = <,91667, 
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Figure 21. 3° x 1,5 km mountain range case. Analytic v profile for 

O = 1,0 (dashed curve) vs. six-layer fine resolution model 
V profile for a = .91667 (solid curve) at 44°N, 
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Figure 22. 3^ x lo5 km mountain range case» Six-layer fine resolution 

model V profiles for O - ,91667 (solid curve) , O - , 58333 
(dashed curve), and a = .25 (dotted curve) at 44^N, 
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Figiire 23. 3® x 1.5 km mountain range case^ Analytic v-profile for 

a = 0 (dashed curve) vs. six-layer fine resolution model 
profile for O = ^08333 (solid curve). Both at 44°N. 
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Figure 24. 3° x 1.5 km mountain range case. Streamlines at a = .75 

from: a. six-layer fine resolution model 
b. two-layer fine resolution model 
(dashed lines indicate the topography) . 
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Figure 25. 3° x 1«5 km mountain range case* v-profiles for O = .75 

at 44°N from six-layer fine resolution model (solid curve) 
and two- layer fine resolution model (dashed curve) , 
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Figure 26. 3® x 1.5 km mountain range case, v profiles for a = ,25 

at 44®N from six-layer fine resolution model (solid curve) 
and two-layer fine resolution model (dashed curve) , 
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5. Air Weather Service 1 

(AWS/TF) 

Scott AFB, Illinois 62225 

6. Captain Harry Hughes 1 

AFIT/CIPF 

Wright-Patterson AFB, Ohio 45433 

7. Lt. Steve Schrobo, USN 1 

1102 Leahy Road 

Monterey, California 93940 

8. Captain John L. Hayes 3 

Air Force Global Weather Central 

PSC #2, Box 7141 

Offutt AFB, Nebraska 68113 
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